Bound-state corrections in laser-induced nonsequential double ionization 
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We perform a systematic analysis of how nonsequential double ionization in intense, linearly 
polarized laser fields is influenced by the initial states in which both electrons are bound, and by 
the residual ionic potential. We assume that the second electron is released by electron-impact 
ionization of the first electron with its parent ion, using an S-Matrix approach. We work within 
the Strong-Field Approximation, and compute differential electron momentum distributions using 
saddle-point methods. Specifically, we consider electrons in Is, 2p, 3p and localized states, which are 
released by either a contact or a Coulomb-type interaction. We also perform an adequate treatment 
of the bound-state singularity which is present in this framework. We show that the momentum 
distributions are very sensitive with respect to spatially extended or localized wave functions, but 
are not considerably influenced by their shape. Furthermore, the modifications performed in order 
to overcome the bound-state singularity do not significantly alter the momentum distributions, 
apart from a minor suppression in the region of small momenta. The only radical changes occur if 
one employs effective form factors, which take into account the presence of the residual ion upon 
rescattering. If the ionic potential is of contact type, it outweighs the spreading caused by a long- 
range electron-electron interaction, or by spatially extended bound states. This leads to momentum 
distributions which exhibit a very good agreement with the existing experiments. 



I. INTRODUCTION 

Within the last few years, nonsequential double ioniza- 
tion (NSDI) in strong, linearly polarized laser fields has 
attracted a great deal of attention, both experimentally 
and theoretically Q]. This interest has been triggered 
by the outcome of experiments in which the momentum 
component parallel to the laser field polarization could be 
resolved, either for the doubly charged ion 0, or for both 
electrons 0- Indeed, the observed features, namely two 
circular regions along the parallel momenta pi\\ = p 2 \\ 

peaked at p x \\ = p 2 \\ — ±2y / £/p, with U p the pondero- 
motive energy, are a clear fingerprint of electron-electron 
correlation, and can be explained by a simple, three-step 
rescattering mechanism , Thereby, an electron leaves 
an atom through tunneling ionization (the "first step"), 
propagates in the continuum, being accelerated by the 
field (the "second step"), and recollides inelastically with 
its parent ion (the "third step" ) . In this collision, it trans- 
fers part of its kinetic energy to a second electron, which 
is then released. 

From the theoretical point of view, there exist mod- 
els, both classical p| | o , Q and quantum-mechanical 
H H [H [H d ESllll , based on such a mechanism, 
which qualitatively reproduce the above-mentioned fea- 
tures. They leave, however, several open questions. A 
very intriguing fact is that, for instance, a very good 
agreement with the experiments is obtained if the inter- 
action through which the first electron is dislodged is of 
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contact type, and if no Coulomb repulsion is taken into 
account in the final electron states. This agreement wors- 
ens if this interaction is modeled in a more refined way, 
considering either a more realistic, Coulomb-type inter- 
action, or final-state electron-electron repulsion. Specifi- 
cally, in recent publications, such effects have been inves- 
tigated in detail using both an S-Matrix computation and 
a classical ensemble model, and have been interpreted in 
terms of phase-space and dynamical effects 0, [f| • This 
analysis has been performed within the Strong- Field Ap- 
proximation (SFA) |15| . which mainly consists in neglect- 
ing the atomic binding potential in the propagation of 
the electron in the continuum, the laser field when the 
electron is bound or at the rescattering, and the internal 
structure of the atom in question. 

Within this framework, the NSDI transition ampli- 
tude is written as a five-dimensional integral, with a 
time-dependent action and comparatively slowly varying 
prefactors. Such an integral is then solved using saddle- 
point methods. Apart from being less demanding than 
evaluating such an integral numerically |3, |9l, o r solv- 
ing the time-dependent Schrodinger equation |l6j. these 
methods provide a clear space-time picture of the phys- 
ical process in question. In particular, the results are 
interpreted in terms of the so-called "quantum orbits". 
Such orbits can be related to the orbits of classical elec- 
trons, and have been extensively used in the context of 
above-threshold ionization, high-order harmonic genera- 
tion h7| and, more recently, nonsequential double ioniza- 
tion mm m [3. 

The fact that, in [T2I Il3j. the crudest approximation 
yields the best agreement with the experiments, seems 
to indicate that the presence of the residual ion, which 
is not taken into account, screens both the long-range 
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interaction which frees the second electron and the final- 
state repulsion. This suggests that the presence of the 
ionic binding potential in the physical steps describing 
nonsequential double ionization, i.e., tunneling, propaga- 
tion and electron-impact ionization, should somehow be 
incorporated. Indeed, in recent studies, it was found that 
Coulomb focusing considerably influences the NSDI yield 

m 

Another possibility is related to how the initial states 
in which the electrons are bound affect the electron- 
momentum distributions. Indeed, the poor agreement 
between the computations with the Coulomb interaction 
and the experiments may be related to the fact that Is 
states have been used in this case, instead of states with 
a different shape or spatial symmetry, such as, for in- 
stance, p states. Furthermore, it may as well be that an 
additional approximation performed in |l2l Il3j for the 
contact interaction, namely to assume that it takes place 
at the origin of the coordinate system, contributes to the 
good agreement between theory and experiments in this 
case. Physically, this means that the spatial extension 
of the wave function of the second electron is neglected. 
Such an approximation has not been performed in the 
computations for the Coulomb interaction discussed in 
0,0], and, up to the present date, there exist no sys- 
tematic studies of its influence in the context of NSDI. 

In this paper, we investigate such effects in the sim- 
plest possible ways. First, we assume that both electrons 
are initially in hydrogenic 2p and in 3p states, instead 
of in s states, as previously done 0, 0, 0, 0. One 
should note that, in contrast to Helium, for which s states 
are more appropriate, p states yield a more realistic de- 
scription of the outer-shell electrons in neon and argon, 
respectively. Since the two latter species are used in most 
experiments, the choice of p states is justified. This is in- 
cluded in the transition amplitude as a form factor, and 
does not modify the saddle-point equations. In both p 
and s - state cases, we consider that the bound-state 
wave function of the second electron is either localized at 
the origin or extends over a finite spatial range, for the 
contact and Coulomb interactions. This provides infor- 
mation on how the initial state of the second electron in- 
fluences electron-impact ionization, and hence the NSDI 
yield. 

A further improvement consists in overcoming the 
bound-state singularity, which is present in the saddle- 
point framework, and which has not been addressed in 
00. 

For this purpose, we use a slightly modified 
action, with respect to that considered in 0, so 
that the tunneling process and the propagation of both 
electrons in the continuum is altered. Such corrections 
depend on the initial wave function of the first electron. 
Hence, they shed some light on how this wave function 
affects the electron-momentum distributions. In particu- 
lar, we investigate how such corrections influence several 
features in the momentum distributions, such as their 
shapes, the cutoff energies or the contributions from dif- 
ferent types of orbits to the yield. 



Finally, we employ a modified form factor for the first 
electron, upon return, which takes into account the ionic 
potential. This is a first step towards incorporating the 
residual ion in our formalism. As it will be discussed 
subsequently, this provides a strong hint that the ion is 
important, in order to achieve a good agreement between 
theory and experiment. 

The manuscript is organized as following. In Sec. [HJ 
we provide the necessary theoretical background for un- 
derstanding the subsequent discussions. In Sees. 11111 II VI 
andlVl we present our results, and, finally, in Sec. I VII we 
state our conclusions. 



II. BACKGROUND 

A. Transition amplitude 

The transition amplitude of the laser-assisted inelastic 
rescattering process responsible for NSDI, in the strong- 
field approximation and in atomic units, is given by 

/OO ft P 

alt / dt' / d 3 kV P]M V K0 eM*S(t,t' ,pj,k)], 
-oo J — oo J 

with the action 



S{t,t', Pj ,k) = - 




~j\k + A(T)} 2 dT+\E Q1 \t' + \E Q2 \t. 

Eq. describes the following physical process: at a 
time t', both electrons are bound (|-Ebi| an d I-E02 1 denote 
the first and second ionization potentials, respectively). 
Then, the first electron leaves the atom by tunneling ion- 
ization, and propagates in the continuum from the time 
t' to the time t, only under the influence of the exter- 
nal laser field E(t) = —dA(t)/dt. At this latter time, 
it returns to its parent with intermediate momentum k, 
and gives part of its kinetic energy to the second elec- 
tron through the interaction V12, so that it is able to 
overcome the second ionization potential | -E7o2 1 - Finally, 
both electrons reach the detector with final momenta Pj 
(j = 1,2). All the influence of the binding potential V 
and of the electron-electron interaction V12 is included in 
the form factors 

Vp J , k =<P2 + A(t),p 1 +A(t)|y 12 |k + A(i),4 2) > (3) 
and 

y k , =<k + A(f')l^l^ 1) >, (4) 

which are explicitly given by 

V K o = 7^372 |d 3 r 1 exp[i(k + A(t')).r 1 ]V(r 1 )4 1) (r 1 ) 

(5) 
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and 



V, 



Pj. k 



d 3 rid 3 r 2 exp[i(pi — k).ri] 



(2tt) 9 / 2 

x exp[i(p 2 + A(t)).r 2 ]V 12 (r 2 , n)^ 2) (r 2 ), (6) 



respectively. The binding potential V(ri) will be taken 
to be of Coulomb type, and the interaction Vi2(r2,ri) 
through which the second electron is released will be cho- 
sen to be of contact or Coulomb type. The initial state 
■00 (ri) of the first electron at the moment of its ion- 
ization will be taken as a hydrogenic s or p state, and 

(2) 

the wave function tp^ (^2) of the second electron at the 
moment of its release is either chosen as a hydrogenic 
state, or a Dirac delta state localized at r 2 — 0. In Eq. 

, we neglect final-state electron-electron repulsion (for 
a discussion of this effect see Refs. [8lll3jh 



B. Saddle-point analysis 

We solve Eq. applying the steepest descent method, 
which is a very good approximation for low enough fre- 
quencies and high enough driving-field intensities. In 
this case, we must find k, t' and t so that S(t, £',p n ,k) 
(n = 1,2) is stationary, i.e., its partial derivatives with 
respect to these parameters vanish. This yields 

[k + A(<')] 2 = -2|£ i| (7) 



J> J+ A(i)] 2 = [k + A(t)] 2 -2|£ 2| (8) 



dr [k + A(r)] = 0. (9) 



Eq. Q gives the energy conservation during tunneling 
ionization, and, for a non-vanishing ionization potential, 
has no real solution. Consequently, t, t' and k are com- 
plex quantities. In the limit |-Eoi| — ¥ 0, Eq. J7J describes 
a classical electron leaving the origin of the coordinate 
system with vanishing drift velocity. Eq. JHJ expresses 
energy conservation at t, in an inelastic rcscattcring pro- 
cess in which the first electron gives part of its kinetic 
energy to the second electron, so that it can overcome 
the second ionization potential and reach the continuum. 
Finally, Eq. Q yields the intermediate electron momen- 
tum constrained by the condition that the first electron 
returns to the site of its release. 

The saddles determined by Eqs. 0-© always occur in 
pairs that nearly coalesce at the boundaries of the energy 
region for which electron-impact ionization is allowed to 
occur, within a classical framework. Such a boundary 
causes the yield to decay exponentially, leading to sharp 
cutoffs in the momentum distributions. 



If written in terms of the momentum components par- 
allel and perpendicular to the laser field polarization, 
Eq. JHJl reads 

2 2 

fell + A(t)] 2 = [k + A(<)] 2 - 2\E a2 \ - Y, Pi± (10) 



and describes a hypersphere in the six-dimensional 
(PjiuPjU.) space. This hypersphere delimits a region in 
momentum space for which electron-impact ionization is 
"classically allowed" , i.e., exhibits a classical counterpart. 
For constant transverse momenta, Eq. IjlUI) corresponds 
to a circle in the (Pin 3P2I1) plane centered at —A(t) and 
whose radius is given by the difference between the ki- 
netic energy £Jki n (i) = 1/2 [k + A(t)] 2 of the first elec- 
tron upon return and the effective ionization potential 
I-E02I = I-E02I + Z) 2 =i P 2 J_/2- Clearly, this radius is most 
extensive if the final transverse momenta Pj±), (J = 1, 2) 
are vanishing, such as in the examples provided in Sec. 
IV. 

In order to compute the transition probabilities, we 
employ a specific uniform saddle-point approximation, 
whose only applicability requirement is that the saddles 
occur in pairs Unless stated otherwise (e.g., in 

Sec. IV), we reduce the problem to two dimensions, using 
Eq. © and the fact that the action J2J is quadratic in k. 
Details about this method, in the context of NSDI, are 
given in [111 1 121 Il4| (for above-threshold ionization 
and high-order harmonic generation, c.f., 0] and [2p| . 
respectively) . 

The momentum distributions of electrons for various 
types of interaction V12 read 



M 



d 2 pi± 



d 2 p 2 ±\M L +M R \ 



(11) 



where the transverse momenta have been integrated over, 
and Ml and Mr give the left and the right peak in the 
momentum distributions, respectively, computed using 
the uniform approximation. We consider a monochro- 
matic, linearly polarized field, so that the vector poten- 
tial reads 



A(t) = — cos(wi)e 2 



(12) 



In this case, M R = M(t, t',p) and M L = M(t - T/2, t' - 
T/2,p), where T = 2ir/ui denotes a period of the driv- 
ing field. We use the symmetry property \M(t, t',p)\ = 
\M(t - T/2, t' - T/2, -p)\ to compute the left peak. One 
should note that, for other types of driving fields, such as 
few-cycle pulses, this condition does not hold and each 
peak must be computed independently 0, Il4| . 



III. INITIAL P STATES 

Within the formalism discussed in the previous section, 
the first and second electron, so far, have been assumed 
to be initially in Is- or zero-range-potential bound states, 
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whose energies |2?oi| and I-E02I are taken to be the first 
and second atomic ionization potential, respectively. In 
most experiments, however, species such as neon and ar- 
gon are used, for which the outer-shell electrons are in 2p 
and 2>p states, respectively. For this reason, such states 
should provide a more realistic modeling of laser-induced 
nonsequential double ionization. For symmetry reasons, 
only the states with magnetic quantum number m = 
will contribute to the yield. 

In this case, the bound-state wave functions of both 
electrons will be given by 



^210 ( r i) = Caio^j exp[-aj-7j-] cos6» 



(13) 



and 



Wio( r j) = CaioJ-j (1 - a. 3 r. 3 /2) exp[-ar,-r,-] cos 9, (14) 



respectively, where a 3 - = y/2\Eaj\(j = 1,2), and C„io 
(where n is the principal quantum number) denote nor- 
malization constants. For comparison, we will also con- 
sider hydrogenic Is wave functions, which read 



(15) 



In Eqs. (|13fl - l|15|) .the binding energies of the first and the 
second electron were chosen as the first and the second 
ionization potentials, respectively. 

The form factors Vp^k, for 2p and 3p initial states, 
read 



V, 



(2p) 



Pj.k 



f?(Pi,k) 



and 



V. 



(3p) 



?7(Pi,k) 



[2|^ 02 

p(p 2 ~2\E 02 \) 



+ ( Pl ^p 2 ) (16) 



[2\E t 



02 



^T + (P1-P2), (17) 



respectively, with p = p x + p 2 — k + A(t). Thereby, 
(Pi P2) means that the momenta of both particles 
are interchanged, and rj(pj,\a) (j = 1,2) is a function 
which depends on the interaction in question. The cor- 
responding form factor obtained for an initial state H15|) 
is given by 



r 



(Is) 



??(Pi,k) 



[2\Eq 



P 2 ] 2 



+ (Pi^P 2 ). (18) 



A. Contact-type interaction 

As a first step, we will assume that the second electron 
is released by a contact-type interaction 



V12 ~ S(ri - r 2 ). 



(19) 



In this case, in Eqs. (|16(l - (|18|l . r)(pj, k) is a constant. The 
differential electron momentum distributions computed 
with such form factors are depicted in Figs. l(a)-l(c), 
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FIG. 1: Electron momentum distributions computed using a 
contact-type interaction, as functions of the electron momen- 
tum components parallel to the laser-field polarization. The 
external field was taken to be monochromatic, with frequency 
lu = 0.057 a.u. and intensity I = 8 X 10 14 W/cm 2 (U p = 1.75 
a.u.). In panels (a), (b) and (c), both electrons are initially 
bound in a Is, 2p, and 3p state [Eqs. 1151 . 11311 and (1141) 1 . 
respectively, whereas, in part (d), the first electron is initially 
in a Is state and the wave function of the second electron is 
localized at T2 = 0. In all situations (even for the 3p - state 
case), the atomic species was taken to be neon (|-Em| = 0.79 
a.u. and | J57o2 1 — 1-51 a.u.), in order to facilitate a clear as- 
sessment of the effects caused by the different initial states. 
The transverse momenta have been integrated over. 



as contour plots in the (pi|i,£>2||) plane. In such com- 
putations, only the pair of orbits for which the electron 
excursion times r = t — t' in the continuum are shortest 
have been employed. As an overall feature, the distribu- 
tions are peaked near piu =P2|| = ^y/Up an d spread in 
the direction perpendicular to the diagonal p x u = p 2 |i- 

An inspection of the form factors l|16 [l -l|18 [l . f° r con " 
stant ?j(pj,k), explains this behavior. Indeed, such form 
factors are large if their denominator is small. Since |-Eq 2 | 
is constant, this condition implies that p = p 1 + p 2 — 
k + A(t) is small. To first approximation, since the first 
electron returns at times t close to the minimum of the 
electric field, one may assume that the vector potential at 
this time and the intermediate electron momentum are 
approximately constant. Furthermore, in the model, the 
field is approximated by a monochromatic wave and k is 
given by Eq. ©. Hence, a rough estimate of these quan- 
tities at the return times yields A(t) ~ 2^JTTp and k ~ 0, 
respectively. Thus, p will be small mainly if pi = — p 2 , 
so that contributions along the anti-diagonal p^ii = — p 2 \\ 
will be enhanced. 
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Such contributions get more localized near the maxima 
for highly excited initial states due to the increase in the 
exponent of the denominator. A direct look at the above- 
stated form factors confirms this interpretation, yielding 
maxima along the anti-diagonal and near ±^/I7p. 

Interestingly, the distributions obtained for the con- 
tact interaction are quite different from the circular dis- 
tributions peaked around pi\\ = p 2 \\ — ±2y / [/p observed 
experimentally. Indeed, in order to obtain such distri- 
butions, it is necessary to assume that the initial wave 
function of the second electron is localized at r 2 = 0. 
This is formally equivalent to taking 

V^r, - r 2 )4 2) (r 2 ) ~ <5(n - r 2 )«5(r 2 ). (20) 

Eq. (|2U[1 yields a constant form factor V Pj ^. In Fig. 

1. (d), we present the distributions computed using Eq. 
H2Up. which exhibit a very good agreement with the ex- 
periments. This means that, in reality, the effective wave 
function of the second electron is very localized, most 
probably due to refocussing 0] , or screening effects (2^ • 

B. Coulomb- type interaction 

We will now consider that the second electron is re- 
leased by a Coulomb type interaction, given by 

^2 = l/|ri-r 2 |. (21) 

In this case, in the form factors H16(l - (|18(1 . T] is given by 

r/(Pj,V = [p ^ k]2 (J = 1,2), (22) 

respectively. This causes the prefactors to be large when 
Pj — k (j = 1, 2) is small, in addition to the case for which 
pCl. The influence of such form factors on the electron 
momentum distributions is shown in Fig. 2. Apart from 
the broadening along p^M = — p 2 u caused by the spatial 
extent of the bound-state wave functions (c.f. Sec. IIII A|) . 
the distributions exhibit maxima near the axispiii = or 
P2\\ = 0. Such maxima are due to the factor (|22() in Eqs. 
(|T6|l - ljT%|l . characteristic of the Coulomb-type interaction, 
which is large for pj ~ k. Since, to first approximation, 
contributions from regions of small k dominate the yield, 
one expects maxima in momentum regions where cither 
Pi or p 2 are small. 

Furthermore, as compared to the yields obtained using 
a Coulomb-type interaction and ls-states, there exists 
a small additional broadening in the distributions, with 
respect to the diagonal pi\\ =p 2 \\,3JS well as an increase in 
the contributions from regions where such momenta are 
small. Such effects get more pronounced as the principal 
quantum number increases, as shown in Figs. 2.(b) and 

2. (c). 

However, such modifications do not alter the distri- 
butions in a significant way. More extreme changes 
occur, for instance, if a contact-type interaction, i.e., 
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FIG. 2: Electron momentum distributions computed using a 
Coulomb-type interaction, as functions of the electron mo- 
mentum components parallel to the laser-field polarization. 
The field and atomic parameters are the same as in Fig. 1. 
In panels (a), (b) and (c), both electrons are taken to be ini- 
tially in a Is, 2p, and 3p state, respectively, whereas in panel 
(d) the first electron is in a Is state, while the spatial ex- 
tension of the bound-state wave function has been neglected. 
The transverse momenta have been integrated over. 

?7(pj,k) = const, in Eq. I|19fl . is taken into account. 
Still, less localized bound states for both electrons will 
cause a broadening in the momentum distributions. If 
the second-electron wave function is localized at the ori- 
gin, the form factor (|18|) reduces to 

Vp " k ~ [ Pl = k p + (Pl ~ P2) - (23) 

The distributions for the latter form factor are displayed 
in Fig. 2.(d). In the figure, one observes a considerable 
reduction of the broadening along the anti-diagonal px\\ = 
—p 2 \\- However, the distributions still exhibit the two 
sets of maxima near the axis pi\\ = or p 2 \\ — 0. This 
is expected, since such maxima are a fingerprint of the 
Coulomb interaction. 

The results in this section show that the shapes of the 
momentum distributions in NSDI are not only influenced 
by the type of interaction by which the second electron is 
dislodged but, additionally, depend on the spatial exten- 
sion of the wave function of the state where it is initially 
bound. In fact, radically different shapes are observed 
if this wave function is either taken to be localized at 
r 2 = or exponentially decaying. This is true both for a 
contact- and a Coulomb-type interaction (c.f. Figs. l.(d) 
and 2.(d)). 

On the other hand, if different initial states are taken, 
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for the same type of interaction, there are no significant 
changes in the shapes of the distributions as long as such 
states extend over a finite spatial range. This is explicitly 
seen by comparing yields obtained using bound states 
with different principal quantum numbers. This is re- 
lated to the fact that the wave functions l|13 |l -l|15 |l were 
chosen such that the bound-state energy always corre- 
sponds to the second ionization potential. Hence, even if 
their shape changes, the spatial extension of such wave 
functions is roughly the same. 

It is still, however, quite puzzling that the best agree- 
ment with the experimental findings occurs for the crud- 
est approximations, both for the interaction and the ini- 
tial bound-state wave function, i.e., for a contact-type 
interaction and a wave function localized at Y2 = 0. 
Indeed, taking either a more realistic type of electron- 
electron interaction, spatially extended bound states, or, 
still, bound states which are, in principle, a more refined 
description of the outer-shell electrons, only worsens the 
agreement between experiment and theory. 

If the main physical mechanism of NSDI is electron- 
impact ionization, there exist two main possibilities for 
explaining this discrepancy. Either the second electron 
is bound in a highly localized state and both electrons 
collide through an effective short-range interaction, as 
the present results suggest, or the tunneling ionization, 
as well as the electron propagation in the continuum, 
must be improved. The first issue may be addressed by 
including the influence of the residual ion in the process, 
whereas the second issue may be dealt with in several 
ways. For instance, in the subsequent section, we will 
consider corrections of a more fundamental nature, which 
alter the semiclassical action and thus the orbits of the 
electrons. 



IV. TREATMENT OF THE BOUND-STATE 
SINGULARITY 

Up to the present section, we have implicitly assumed 
that the form factors V Vj ^_ and Vk,o are free of singu- 
larities and slowly varying in comparison to the time- 
dependent action. However, this is not always true. In- 
deed, in the saddle-point framework, the form factor Vk,o 
is singular if the electron is initially in a state described 
by an exponentially decaying wave function, such as Eqs. 
(fH?|l - <(T5jl . More specifically, in this case, 



Vk,o oc 



f(k + A(t')) 



[k + A(t')f + 2\E Q 



(24) 



where n is an integer number. In this case, according to 
Eq. (JZJ), the denominator vanishes. Due to this singular- 
ity, this form factor does not vary slowly with respect to 
the semi-classical action and thus must be incorpo- 
rated in the exponent. Therefore, we take the modified 
action 

S(t,t', Pi ,k) =S(i,i / ) p i ,k)-iln[y k ,o] (25) 



in the transition amplitude Q. This causes a change in 
the first and third saddle-point equations, which will de- 
pend on the initial bound state in question. In particular, 
we will consider that the first electron is initially in the 
hydrogenic states Is, 2p, and 3p. This is a legitimate as- 
sumption, since the binding potential of a neutral atom, 
from which the first electron tunnels out, is of long-range 
type. For the states Is, 2p, and 3p, Vk^o reads 



V, 



(Is) 



V2 (2|£ i|) 5/2 



k.O 



2|E 



01 



V, 



(2p) 



k.O 



7T (v2+2|£ l|) 2 ' 



and 



V, 



( 3P ) _ 8*( v /2l^) 5 / 2 i;(v 2 -2|£ 01 |) 



k.O 



V3~7r(v 2 +2|£oi| 



(26) 



(27) 



(28) 



respectively, where v = [k + A(i')] denotes the initial 
electron drift velocity. The explicit expressions for the 
saddle point equations then become 



[k + A(f)] = -2|£ i| + C(*. O-E(i') (29) 



and 



J dr[k + A(T)]+C(fc,t') = 0, 



(30) 



respectively, where £(fc,i') = — idk ln[Vt,o] is a correction 
which depends on the initial bound state. Thus, there is 
an effective shift in the ionization potential at the tun- 
neling times, and a modification in the return condition. 
Consequently, the orbits change. Apart from that, from 
the technical point of view, the transition amplitude is 
no longer reducible to a two-dimensional integral, so that 
the problem is far more cumbersome. 

The modifications in the equation describing tunnel- 
ing ionization allow the existence of solutions for which 
Rc[v] 7^ 0. This did not occur in Eq. 0, for which this 
quantity was purely imaginary, and, physically, means 
that there are in principle changes, maybe even enhance- 
ments, in the probability that the first electron tunnels 
out at t'. 

Furthermore, Eq. I|3U|I . if written in terms of the com- 
ponents of the intermediate momentum k parallel and 
perpendicular to the laser field polarization, has, apart 
from the trivial solution kj^ = 0, additional solutions 
for which k^ =/= 0. Thus, in principle, the first electron 
may have, during the tunnel ionization and upon return, 
a non-vanishing drift velocity component transverse to 
the laser-field polarization. We regard this possibility, 
however, as non-physical, and therefore will mainly con- 
centrate on the case of vanishing k^. Despite of that, the 
results obtained for nonvanishing k^will be briefly dis- 
cussed in Sec. IIVBI For the return condition (JSJ) this is 
not possible and k^ is always vanishing. In the follow- 
ing, we will investigate how the corrections in the action 
affect the momentum distributions. 



7 







A. Vanishing kx 

In this section, we will consider that the first elec- 
tron has vanishing intermediate momentum components 
k±. Physically, this means that the dynamics of NSDI 
is mainly taking place along the laser field polarization, 
which is the intuitively expected situation. In Fig. 3, we 
present the electron momentum distributions computed 
employing the modified saddle-point equations and the 
action 1)25(1 ■ for the same initial states and types of as 
in Figs. 1 and a contact-type interaction. In general, 
the distributions in Fig. 3 are very similar to the former 
ones, with, however, a suppression in the region of small 
parallel momenta. This is true even if different correc- 
tions are taken into account, as it is the case if the first 
electron is initially in a Is, 2p and 3p state (Figs. 3. (a), 

3. (b) and 3.(c), respectively). In the specific case of a 
localized bound-state wave function for the second elec- 
tron, there is also a minor displacement of the maxima 
towards smaller parallel momenta(c.f. Fig. 3(d)). 

The suppression persists if the second electron is re- 
leased by a Coulomb-type interaction, as shown in Fig. 

4. Specifically for this interaction, the corrections lead 
to a suppression of the secondary maxima in the small- 
momentum region, which were present in Fig. 2. 

In the following, we will analyze these differences in 
terms of the so-called quantum orbits, obtained by solv- 
ing the saddle point equations. We will consider both 
the saddle-point equations in the presence and absence 
of corrections to the bound-state singularity, i.e., Eqs. 
Pty. © and (J3DJ, and 0-©, respectively. We re- 
strict ourselves to vanishing final transverse momenta 
and longitudinal momentum components along the diag- 
onal pi I = Pin = P2\ | ■ F° r this particular case, the energy 
region for which electron-impact ionization is classically 
allowed is most extensive. 

In Fig. 5, we display the solutions of the saddle-point 
equations for the rescattering times t and the interme- 
diate momentum k. The upper and lower panels in the 
figure give the real and imaginary parts of such variables, 
respectively. The real parts of t and k correspond to the 
solutions of the equations of motion of a classical elec- 
tron in an external laser field, and almost merge at two 
distinct parallel momenta. These momenta are related 
to the maximal and minimal energy for which the second 
electron is able to overcome \Eq2 | . Beyond such momenta, 
there are cutoffs in the distributions, and the yield decays 
exponentially. The imaginary parts of such variables are 
in a sense a measure of a particular physical process be- 
ing classically allowed or forbidden. Indeed, the fact that 
|Im[i]| and |Im[fc]| are vanishingly small between the min- 
imal and maximal allowed momenta are a consequence 
of both electron-impact ionization and the return con- 
dition being classically allowed in this region. As the 
boundaries of this region are reached, |Im[t]| and |Im[fc]| 
increase exponentially. Interestingly, both the real and 
imaginary parts of such variables, as well as the cutoff 
momenta, remain practically inaltered upon the changes 
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FIG. 3: Electron momentum distributions computed with us- 
ing contact-type interaction, as functions of the electron mo- 
mentum components parallel to the laser-field polarization, 
for the same field and atomic parameters as in the previous 
figures. We introduce corrections to the bound-state singu- 
larity by employing the modified action l|25|l and saddle-point 
equations, taking the solutions displayed in Figs. 6. (a) and 
6.(c). In panels (a), (b) and (c), both electrons are taken to 
be initially in a Is, in a 2p, and in a 3p state, respectively, 
whereas in panel (d) the first electron is initially in a Is state, 
while the spatial extension of the bound-state wave function 
of the second electron has been neglected. The transverse 
momenta have been integrated over. 



introduced in this section. This is not obvious, since the 
bound-state corrections in question alter the return con- 
dition [c.f. Eq.lHU]. 

There exist, however, modifications in the tunneling 
times t' , which are explicitly shown in Fig. 6. Specif- 
ically, the corrections in the tunneling condition, which 
leads to Eq. I|29() . cause a splitting in the solutions of 
Eq. 0. This follows from the fact that small variations 
in the stationary-action trajectories contributes quadrat- 
ically to S(t, i', pj, k) and Vk,o, so that S(t, £',Pj,k) at- 
tains two stationary trajectories for each of the former 
ones. Strictly speaking, a similar splitting also occurs 
for k and t. In practice, however, the difference between 
the two different sets of solutions is vanishingly small, 
and thus not noticeable in Fig. 5. The different sets 
of solutions are depicted in Figs. 6(a) and 6.(c), and 
6.(b) and 6(d), respectively. The real parts Re[i'] ex- 
hibit only minor differences, with occur for the shorter 
orbits and small momenta and eventually disappear as 
the upper cutoff is approached. Depending on the type 
of correction, such times either distance themselves from, 
or become slightly closer to the peak-field times (Figs. 
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FIG. 4: Electron momentum distributions computed using 
the Coulomb-type interaction, as functions of the electron 
momentum components parallel to the laser-field polariza- 
tion, for same field and atomic parameters as in the previous 
figure. We introduce corrections to the bound-state singu- 
larity by employing the modified action l|25|l and saddle-point 
equations, and the solutions in Figs. 6(a) and 6.(c). In panels 
(a), (b) and (c), both electrons are taken to be initially in a 
Is, in a 2p, and in a 3p state, respectively, whereas in panel 
(d) the first electron is initially in a Is state, while the spa- 
tial extension of the bound-state wave function of the second 
electron has been neglected. The transverse momenta have 
been integrated over. 



6. (a) and 6.(b), respectively). Thus, one could expect an 
enhancement in the contributions from the shorter orbits 
near the origin of the (pi|i,P2||) plane, in the former case, 
and a suppression in the latter case. However, we have 
used the solutions in Fig. 6.(b) and (d) for computing 
the contour plots in Figs. 3 and 4, and obtained a sup- 
pression in the yield. This is a clear indication that the 
changes in Im[i'] and in the time-dependent action play 
a more important role than those in Re[i']. 

In Figs. 6.(c) and 6.(d), we present the imaginary parts 
of t' , which clearly shift towards smaller, and larger val- 
ues, respectively, when the corrections <,(k,t') are taken 
into account. The higher the initial state lies, the larger 
such shifts are. Physically, there exists a correspondence 
between such imaginary parts and the probability that 
the first electron tunnels out and reaches the contin- 
uum. This means that, by using a slightly modified ac- 
tion in order to overcome the Coulomb singularity, one 
is changing the effective potential barrier at t' for the 
first electron. In general, such a barrier has a significant 
influence on the distributions. Indeed, recently, we have 
shown, within the context of nonsequential double ion- 
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FIG. 5: Rescattering times, together with the intermediate 
momentum k\\, as functions of the parallel momentum p|| 
along the diagonal p^i = p 2 ||- The final transverse momenta 
Pj±(j = 1, 2) and the intermediate transverse momentum kx 
are taken to be vanishing. The real and imaginary parts of 
such quantities are displayed in the upper [(a), (b) and (c)] 
and lower [(d), (e) and (f)] panels, respectively. The field 
and atomic parameters are the same as in the previous fig- 
ure. The uncorrected variables are given by the thick light 
gray curves in the figure, while the variables with corrections 
corresponding to initial Is, 2p and 3p states are given by the 
gray, dark gray and black curves, respectively. The longer 
and the shorter orbit are indicated by dashed and solid lines, 
respectively. 



ization with few-cycle laser pulses, that the importance 
of the contributions of a particular orbit or set of orbits 
to the yield is highly dependent on |Im[t']|. The smaller 
this quantity is, the larger is the tunneling probability for 
the first electron [l4|. As a direct consequence, contri- 
butions from orbits with small |Im[t']|, i.e., with a large 
tunneling probability, dominate the yield. In the present 
case, however, since both orbits are being equally shifted, 
this should not influence the distributions qualitatively. 
One should note that, even in the momentum region for 
which electron-impact ionization is allowed, |Im[i']| is al- 
ways nonvanishing. This is a direct consequence of the 
fact that tunneling ionization is a classically forbidden 
process. 

Subsequently, we compute the counterparts of Fig. 3 
and 4 (Figs. 7 and 8) using the solutions displayed in 
Fig. 6.(b) and 6.(d). Also in this case, in general, there 
is a suppression in the yield in the region of small parallel 
momenta, with, however, a slightly different substructure 
in the Coulomb-interaction case. 



B. Nonvanishing kx 

The modifications introduced in the return condition 
for the first electron [Eq. (|30|l ] allow the intermediate 
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FIG. 6: Tunneling times as functions of the parallel momen- 
tum p|| along the diagonal = p 2 ||- The final transverse 
momenta Pj±(j = 1,2) and the intermediate transverse mo- 
mentum kx are taken to be vanishing. The real and imagi- 
nary parts of such quantities are displayed in the upper [(a) 
and (b)] and lower [(c) and (d)] panels, and the same pairs of 
orbits are depicted in the left [(a) and (c)] and right [(b) and 
(d)] panels, respectively. The field and atomic parameters are 
the same as in the previous figure. The uncorrected variables 
are given by the thick light gray curves in the figure, while the 
variables with corrections corresponding to initial Is, 2p and 
3p states are given by the gray, dark gray and black curves, 
respectively. The longer and the shorter orbit are indicated 
by dashed and solid lines, respectively. 



momentum k to have a nonvanishing component per- 
pendicular to the laser-field polarization. This implies 
that the first electron, during tunneling ionization and 
when it returns, is being deviated from its original di- 
rection. Although such an effect is unphysical, we will 
briefly discuss its consequences. For that purpose, we 
will consider the simplest corrections to the bound-state 
singularity discussed in this paper, namely those for Is 
initial states. If Eq. (|3U[) is written in terms of the 
intermediate-momentum components k^ and k\ \ perpen- 
dicular and parallel to the laser-field polarization, this 
equation reads 



k\\(t-t')- I A(s)ds+ 



and 



2i[k n +A(t')] 
[2|£ 01 |+ki + [fc||+A(£')] 2 ] 

(31) 



t' + 



2i 



[2|£ 01 |+ki + [fc||+A(*')F 



0, (32) 



respectively. Apart from the trivial solution k^ = 0, the 
condition \.\'2\) can be satisfied by nonvanishing values of 
this variable. One should note that, in the case with- 
out corrections, this does not hold and only the trivial 
solution exists. 

Fig. 9 depicts the tunneling and rescattering times 
for this case, together with the perpendicular and par- 
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FIG. 7: Electron momentum distributions computed with the 
contact-type interaction, as functions of the electron momen- 
tum components parallel to the laser-field polarization, for 
the same field and atomic parameters as in the previous fig- 
ure. We introduce corrections to the bound-state singularity 
by employing the modified action l|25|l and saddle-point equa- 
tions, taking the solutions in Figs. 6.(b) and 6.(d). In panels 
(a), (b) and (c), both electrons are taken to be initially in a Is, 
in a 2p, and in a 3p state, respectively, whereas in panel (d) 
the first electron is in a Is state, and spatial extension of the 
bound-state wave function of the second electron has been ne- 
glected. The transverse momenta have been integrated over. 



allel components of k. The real parts of such variables 
correspond, as in the previous cases, to a longer and a 
shorter orbit. The momenta, however, for which such 
orbits nearly coalesce, are radically different from those 
in the previous cases discussed in this paper. This is 
due to the fact that a nonvanishing k^ also affects the 
rescattering condition ijHJ), which now reads 



£ | + A(t)] 2 = [k\\ + A(t)] 2 +kl-2\E a2 \-J2p 



For constant final transverse momenta Pj±(j — 1,2), this 
equation describes a circle centered at —A(t) whose ra- 
dius has been altered in k^_. Since, as shown in Figs 9(a)- 
9(d), this radius decreased, kj^ is expected to be almost 
purely imaginary. This is indeed the case, as can be seen 
comparing panels (d) and (h) in the figure. The imag- 
inary parts of such variables also behave following the 
same pattern as previously, growing vary rapidly at the 
momenta for which the real parts approach each other, 
and remaining nearly constant in-between. Interestingly, 
lxsx[t ] is vanishing in this region. This feature is in clear 
contradiction with the fact that tunneling is a process 
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FIG. 8: Electron momentum distributions computed using 
the Coulomb- type interaction, as functions of the electron mo- 
mentum components parallel to the laser-field polarization, 
for same field and atomic parameters as in the previous fig- 
ure. We introduce corrections to the bound-state singularity 
by employing the modified action 125H and saddle-point equa- 
tions, and the solutions in Figs. 6.(b) and 6.(d). In panels 
(a), (b) and (c), both electrons are taken to be initially in a 
Is, in a 2p, and in a 3p state, respectively, whereas in panel 
(d) the first electron is initially in a Is state, while the spa- 
tial extension of the bound-state wave function of the second 
electron has been neglected. The transverse momenta have 
been integrated over. 



which is always forbidden, and therefore requires a non- 
vanishing lm[t'] (c.f. Fig. 5.(d)). For this reason, we will 
not use solutions with nonvanishing for computing 
electron momentum distributions. 



FIG. 9: Tunneling [Panels (a) and (e)] and rescattering times 
[Panels (b) and (f)], together with the parallel [Panels (c) and 

(g) ] and perpendicular components of the intermediate mo- 
mentum k [Panels (d) and (h)], as functions of the parallel 
momentum pn along the diagonal put = p 2 n. The final trans- 
verse momenta Pj±(j = 1, 2) are taken to be vanishing. The 
real and imaginary parts of such quantities are displayed in 
the upper [(a), (b), (c) and (d)] and lower [(e), (f), (g) and 

(h) ] panels, respectively. The field and atomic parameters 
are the same as in the previous figure. The corrected and 
uncorrected yields are given by the black and gray curves in 
the figure, respectively. The tunneling and rescattering times 
are multiplied by u and the intermediate momenta divided 
by y/ui, respectively, so that a direct comparison with Figs. 5 
and 6 can be performed. 



In this context, both the effective charge and a contact- 
type interaction are justified by the fact that the remain- 
ing electrons are screening the charge and the long-range 
tail of the binding potential. 

In Eq. the form factors are given by 



2r?(pi,k) 



[2|£ 02 | + (p 2 + A(i)) 



212 



Pi 



P2), 

(35) 



V. INFLUENCE OF THE ION 

In this section, we take a first step towards including 
the residual ion in our formalism. For that purpose, we 
consider an effective interaction V — V%2 + V- lon at the 
time the first electron returns, where V lon is the ionic 
potential. Physically, this means that the first electron 
interacts not only with the electron it releases, but, ad- 
ditionally, with the residual ion. We take this potential 
to be of either Coulomb or contact type, and assume 
that only the two active electrons contribute to the ionic 
charge. Thus, explicitly, Vi on reads 



V, 



(C) _ 



or 



V, 



-2/N 



-25(ri). 



(33) 



(34) 



(2p) 



27 ? (p 1 ,k) V /(P2+A(t))2 



(pa+A(i)) a P 



+ (Pi <-> P2), 
(36) 

for Is and 2p states, respectively. The prefactor V12 is of 
contact or Coulomb type. Furthermore, for a Coulomb 
or contact ionic potential, 77 is either constant or given 
by Eq. (|22|l . respectively. In order to simplify the com- 
putations, and since only minor differences have been 
observed in this case, we use the model in Sec. II, in- 
stead of the more rigorous approach of Sec. IV in the 
subsequent figures. 

Fig. 10 depicts how the ion affects the electron mo- 
mentum distributions, if its potential is assumed to be 
of Coulomb form (Eq.(J23J)- The upper and lower pan- 
els have been computed for V12 of contact and Coulomb 
type, respectively. In the figure, the distributions resem- 
ble those obtained for the Coulomb-type interaction, if 
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FIG. 10: Electron momentum distributions computed with 
the modified form factors (1351 . 1136^ . as functions of the elec- 
tron momentum components parallel to the laser-field polar- 
ization, for same field and atomic parameters as in the previ- 
ous figures. We consider that the ionic potential is of Coulomb 
type (Eq. |33| V and use the saddle-point model without mod- 
ifications (c.f. Sec. II). In the upper and lower panels, the 
electron-electron interaction V12 was assumed to be of contact 
and Coulomb type, respectively. In panels (a), (b), both elec- 
trons are taken to be initially in Is states, whereas in panels 
(c) and (d), they are initially in 2p states. The transverse 
momenta have been integrated over. 
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FIG. 11: Electron momentum distributions computed with 
the modified form factors I35L136H . as functions of the elec- 
tron momentum components parallel to the laser-field polar- 
ization, for same field and atomic parameters as in the previ- 
ous figures. We consider that the ionic potential is of contact 
type (Eq. (|34|> ^ . and use the saddle-point model without mod- 
ifications (c.f. Sec. II). In the upper and lower panels, the 
electron-electron interaction V12 was assumed to be of contact 
and Coulomb type, respectively. In panels (a), (b), both elec- 
trons are taken to be initially in Is states, whereas in panels 
(c) and (d), they are initially in 2p states. The transverse 
momenta have been integrated over. 



the second state is in a localized state (Figs. 2.(d), 4.(d), 
and 8.(d)). This holds both for Is and 2p initial elec- 
tron states. An inspection of Eqs. 1)35(1 and (|36[1 explains 
this shape. Indeed, in both equations, the functional 
form of T)(j>j, k), which is characteristic of long-range in- 
teractions, favors unequal momenta, leading to patterns 
similar to those observed in Figs. 2, 4, and 8. Further- 
more, in the second terms in V^k, the denominators are 
small if Pj ~ — A(t). Thus, since, A(t) ~ 2y/Up, we 
expect the form factors (|35fl and to be large near 
Vi\\ = V\\\ — ±2^/T7p. Consequently, the yield in the di- 
agonal gets enhanced. This is a feature shared with the 
limit for localized wave functions, so that the distribu- 
tions are similar. 

The subsequent figure (Fig. 11) is the counterpart of 
Fig. 10 for a contact-type ionic potential (Eq. (|34(l ). In 
this case, for all types of electron-electron interaction V12 
and initial bound states, the distributions are strongly 
localized near p 2 || = V\\\ = ±%\/Upi even though their 
shapes are slightly different. This happens due to the fact 
that, in this case, r](pj,]s.) = const. Therefore, the sec- 
ond terms in the form factors (|35|) and i|3o|) are large near 
P2\\ = Pr|| — i2y / C/p, but, in contrast to the Coulomb- 



potential case, no unequal momenta are favored. This 
means that the inclusion of a short-range ionic poten- 
tial leads to a radical improvement in the agreement be- 
tween theory and experiment. In this context, for both 
the contact- and Coulomb- type interactions V12, circu- 
lar shapes reminiscent of those in the experiments are 
only obtained if we consider 2p states. As previously dis- 
cussed, such states provide a more realistic description 
of the outer-shell electrons in neon, as compared to Is 
states. 



VI. CONCLUSIONS 

In this article, we have introduced several technical 
modifications in an S-Matrix theory of laser-induced non- 
sequential double ionization (NSDI), within the Strong- 
Field Approximation, in which this phenomenon is mod- 
eled as the inelastic collision of an electron with its parent 
ion. Such modifications include different initial bound 
states for the first and second electron, an adequate treat- 
ment of the bound-state singularity which exists in our 
framework, and an effective form factor which incorpo- 
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rates the residual ion. We performed a systematic analy- 
sis of their influence on the differential electron momen- 
tum distributions as functions of the parallel-momentum 
components Pju, (j — 1,2) of both electrons. 

Specifically, we consider that the second electron is dis- 
lodged by a contact- and Coulomb-type interaction, and 
assume that both electrons are initially in a Is, 2p or 3p 
hydrogenic state. As an additional case, we assume that 
the first electron is initially bound in a Is state, and that 
the initial wave function of the second electron is local- 
ized at r 2 = 0. For the first electron, we take into account 
only the hydrogenic states, since a neutral atom, in con- 
trast to a singly ionized atom, has a long-range binding 
potential. 

Concerning the initial bound-state wave function of 
the second electron, our results show that the NSD1 
momentum distributions are very sensitive to its spa- 
tial extension, but not to its shape. Indeed, a spa- 
tially extended wave function causes a broadening in the 
electron momentum distributions along the anti-diagonal 
Pi 1 1 — ~~ P2||i even if the second electron is dislodged by a 
contact-type interaction. Circular-shaped distributions, 
as reported in [T2I and observed in experiments Q, Q 
are only obtained for a contact-type interaction under 
the additional condition that the bound-state wave func- 
tion is localized at the origin of the coordinate system, 
i.e., at T2 — 0. In addition to this broadening, if the sec- 
ond electron is released by a Coulomb-type interaction, 
there is an enhancement in the contributions near the 
axis =0 or p 2 \\ — 0. 

All the distributions investigated in this article, how- 
ever, change in a less radical fashion if the second electron 
is taken to be in a a Is, 2p or 3p hydrogenic state, as long 
as they exhibit a spatial extension. In fact, although spe- 
cific changes are observed, such as an additional substruc- 
ture for a Coulomb-type interaction, or more localized 
distributions for a contact-type interaction, the overall 
shapes of such distributions remains similar. 

Furthermore, if the form factor o, which, within our 
model, contains all the influence of the initial state of 
the first electron, is incorporated in the time-dependent 
action, the only noticeable effect is a suppression in the 
yield, for regions of small parallel momenta. Indeed, the 
distributions retain their shapes even if the saddle-point 
equations are modified in this way. Such changes have 
been introduced in order to correct a singularity which 
exists for such the prefactor Vk t o, within the saddle-point 
framework, if the initial bound state is exponentially de- 
caying. 

Finally, the inclusion of the ionic potential at the time 
of rescattering, as the modified form factors V Pjt u, sheds 
some light on why, in the absence of the ion, a contact- 
type interaction localized at the origin of the coordinate 
system yields the best agreement with the experimental 
findings. 

In fact, the ionic interaction leads to form factors which 
are very large near pi\\ = p 2 \\ = ±y^Vp. This causes an 
enhancement in the distributions in this region. If the 



ionic potential is of Coulomb type, this effect is over- 
shadowed by the fact that n(pj,]s.), given by Eq. 
favors unequal momenta. By contrast, if the ionic po- 
tential is given by Eq. (|34f> . which is a good approxi- 
mation for a short-range interaction, n(pj,\i) = const. 
and the enhancement at the diagonal prevails. On the 
other hand, in Sec. Ill and IV, if Vp . ^ = const, (i.e., for 
V12 of contact-type and a localized state for the second 
electron), the very same effect is caused by integrating 
over the phase space. Interestingly, if, in the presence 
of the ion, we consider 2p states, which are more realis- 
tic assumptions for our model, the agreement with the 
experimental findings improves even more. 

In conclusion, the present results indicate that the 
ionic potential is an important ingredient for a realistic 
modeling of NSDI. Indeed, of all the technical modifica- 
tions considered in this paper, which aimed at making 
the model more realistic, this was the only which played 
a major role in improving the agreement between theory 
and experiment. The other modifications either wors- 
ened this agreement, or had almost no influence on the 
momentum distributions. This supports the hypotheses 
raised in previous studies 0, , that the residual ion 
might be screening both the long-range of the Coulomb 
interaction, or the final-state Coulomb repulsion, so that, 
effectively, the electron-electron interaction is of contact- 
type, and the bound-state wave functions are localized. 

We would like to stress out, however, that the treat- 
ment performed in Sec. V is only a first approximation 
for a rigorous study of the ionic potential. There exist, 
in principle, more rigorous methods for incorporating the 
residual ion. The first approach would be to consider the 
ion as a further interaction in our model, and modify 
the transition amplitude accordingly. This is, however, a 
highly non-trivial task, since it would lead to one more 
rescattering and a further integral in the transition am- 
plitude. Another possibility would be to incorporate the 
ionic potential in the propagation of both electrons in 
the continuum. This would allow a clear assessment of 
the Coulomb focusing, which, again, owes its existence to 
the presence of the ion. Indeed, this effect may as well be 
compensating the broadening caused by initial spatially 
extended wave functions. Definite statements on this is- 
sue, however, require a theoretical approach beyond the 
Strong-Field Approximation. 
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